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, We reconsider the issue of Berezinskii-Kosterlitz-Thouless (BKT) transition into an insulating 

' state in the Coulomb-dominated Josephson junction arrays. We show that previously predicted 

picture of the Cooper-pair BKT transtion at T = T2 is valid only under the condition that T2 is 
considerably below the parity-effect temperature T* ~ 0.1 A and even in this case it is not a rigorous 
phase transition but only a crossover, whereas the real phase transition takes place at T\ ~ T2/4. 
Our theory is in agreement with available experimental data on Coulomb-dominated Josephson 
arrays and also sheds some light on the origin of unusual reentrant temperature dependence of 
OO ' resistivity in the array with nearly-criticial ratio Ec / Ej . 



1. Two-dimensional arrays of micron-scale superconducting islands are extensively studied during last years, both 
experimentally and theoretically |5|-Q] ■ It is well-established now that their low-temperature behaviour is deter- 
mined by the competition between Josephson coupling energy Ej and effective charging energy Eq — e 2 /2C, where 
C is some relevant electric capacitance (to be discussed below) . Macroscopic superconductive coherence was observed 
at low tempratures in the arrays with Ej 3> Ec, whereas arrays with Ej <C Ec show insulating behaviour at T — > 0. 
At nearly-critical value of the ratio x = Ej/Ec ~ x cr direct transition between superconductive (SC) and insulating 
(I) behavior as function of x was observed in zero magnetic field [|). MOrevoer, very weak magnetic field B < 1G 
(producing small fractions of flux quantum per unit cell of the array) was shown to switch arrays with x w x cr between 
SC-like and I -like behavior as function of temperature; recently very interesting intermediate region was found 
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where resistance R(T) is basically constant in the temperature range WmK < T < 200mK, which indicates the 
existence of a "2D metal" state sandwiched between SC and I phases. 

The above-mentioned basic properties of 2D arrays are in qualitative agreement with available theories J5|,0 (except 
for the recently observed 2D metal state); however several important features are not understood yet. In particular, 
resistance of the insulating arrays shows purely activated behaviour R(T) oc exp(E a /T) with constant activation 
energy E a through the whole temperature interval studied , whereas theoretically the charge binding Berezinskii- 
t^J- ■ Kosterlitz-Thouless (BKT) transition H|| from the conducting to insulating phase is expected to occur || at the 
C \ temperature T2 ~ Ec/ft- Such a transition should occur due to nearly-logarithmic form of Coulomb interaction 
• between charges in the arrays with self-capacitance of islands Co very small compared to the inter-island (junction) 
capacitance C. In the currently studied arrays the ratio C/Co ~ 100 (as measured at very low temperatures, about 
WmK, cf. e.g. Q), which should result in the logarithmic interaction throughout the whole array (effecitve length of 
I interaction A should be estimated with the account of 3D nature of electric field, which leads P|lC|] to A ~ C/C'o ~ 100) 
■ and, consequently, to the charge binding BKT transition. 

In the case of islands in SC state and under the condition Ej <C Ec the temperature of this transtion was estimated 
P| as T 2 w Ec/'K, whereas in the case of normal islands (i.e. with superconductivity suppressed by magnetic field) it is 
expected to be 4 times lower: T\ « Ec/^K due to the twofold decrease of an elementary charge available. Nevertheless 
no indication of such a transtion in array of SC islands was found experimentally (except in very recent preprint [fj"l|| , 
which is discussed below). Another surprising feature observed in |g{ is nonmonotoneous ("reentrant") temperature 
O ■ behaviour of resistance R(T) of the array with a nearly-critical Ej /Ec ratio at T < 200mK. 

In the present Letter we show that the above experimental observations can be naturally understood once the 
temperature dependence of effective Coulomb interation between the charges in the array is taken into account 
properly. In a very broad sense our analysis follows the ideas of Efetov []12[ who established the background for 
the description of quantum fluctuations in granular superconductors; namely, we consider screening of Cooper-pair 
Coulomb interaction by normal quasiparticles existing in each superconductive island at finite temperatures. However 
we believe that Efetov's treatment of the effect he proposed technically was not quite correct, thus we present here 
another theoretical approach to the same problem. 

Our main qualitative result can be formulated as follows: at the temperatures above the so-called parity-effect 
pMl temperature T* « A/ In M < A [where M = Vv(0)V&nTA - 10 4 - 10 5 , V is the volume of the island and 



i^(0) is the density of states at the Fermi level in absense of superconductivity] the presence of thermal quasiparticles 
[with the number ~ M. exp(— A/T) ^> 1] on each island excludes any possibility to observe Cooper-pairs BKT transtion 
at T2. Since in the most of arrays studied till now the above-defined characterisitic temperature T2 was in the range 
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0.3-0.5 K, whereas parity effect temperature T* w 0.2K < T2, the absense of anything like BKT transtion near T2 is 
quite natural (measurements below T* were not possible in these arrays [||^] since R(T) becomes immeasurably high 
[> 10 9 Ohm]). On the other hand single-electron BKT transtion is a completely different issue: we do expect such 
a transition to be observable at approximately the same temperature T% ~ T2/4 as in the arrays with islands in the 
normal state. 

2. We proceed now to the derivation of our results, and will follow, with one important modification, Ref. ||. In 
the limit when charge tunneling is weak and is important only for the establishment of thermodynamic equilibrium 
an array of superconducting islands can be described by a classical partition function of a form: 
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The first term in the exponent in Eq. (|l|) stands for the electrostatic energy of the array which in the case when only 
mutual capacitance of nearest islands C is of importance corresponds to the logarithmic interaction of the charges in 
two-dimensional array: 
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whereas the second term describes the dependence of a free energy of a superconducting island on the parity of the 
number of electrons n, on this island [113 15 lql. 



The free energy difference D(T) between the islands with odd and even number of electrons can be expressed as 

D{T) = -Tmtanh(f2 oe /T) ( 3 ) 

where Q oe = —T\n(Z dd/Z) and Z Q dd is the "odd grand canonical partition function" introduced in Ref. |l3| for the 
study of parity effect. Accurate expression for the function Q oe (T) can be found in jl5), but for T C A a good 
approximation is given by 

n oe (T)/T sa Me- A/T ; M = Vi/(0)V8ttTA (4) 

The ratio n oe /T in that limit is proportional to the number of thermally exited quasiparticles on one island. 

In terms of statistical mechanics partition function (Q) defines a lattice Coulomb gas in which the fugacities of odd 
charges Y= exp(— D/T) differ from the fugacities of even charges (which are equal to one). Comparison of Eq. ^ with 
Eq. (||) shows then that for T«T* = A/ \nM the parity-dependent free energy difference D(T) sa A - Tin M > T 
and 7< 1, whereas in the opposite limit T ^> T* the quantity D (T) becomes exponentially small and Y is very close 
to one. 

Previously it has been assumed |5| that in the regime when island charges behave as classical variables the main 
difference between the array of normal islands and the array of superconducting islands is that in the array of normal 
islands the charge of each island is quantized in units of e, whereas in the array of superconducting islands the charge is 
quantized in units of 2e. The consequence for the array the electrostatic properties of which are dominated by mutual 
capacitance of nearest neighbours is that the temperature T2 of the BKT transition in the array of superconducting 
islands (appearance of free double charges) should be exactly four times higher than the temperature T\ ~ Ecl^n of 
the BKT transition in the analogous array of normal islands (appearance of free single charges). Comparison with 
Eq. (Q) shows that such description of the array of superconducting islands would be correct only in the limit of 
D(T)/T — > 00. Since D(T) is always finite this description turns out to be misleading. The behaviour of the array at 
temperatures close to T2 = Ec /ft depends qualitatively on the relation between T2 and T* ; we consider both cases in 
turn. 

3. At T2 > T* an array of the superconducting islands is described by practically the same partition function as an 
array of the normal islands, since DiTz) « T2 in that case. The phase transition into insulating state in such system 
can be associated with the binding of the charges ±1 into neutral pairs; it takes place at the temperature T\ which 
is slightly lower than the simple estimate — Ec/4:ir which can be obtained by comparison of the single charge 
energy with its entropy. The difference between Ti and is related to the renormalization of charge interaction by 
bound pairs of charges and decreases with decrease in fugacities. The appearance of the free single charges induces 
the screening of the Coulomb interaction for all types of charges and therefore the double charges also are free at 
T > Ti. Not even a trace of a separate phase transition related to debounding of double charges can be expected to 
be observed in such a situation, which was realized in the experiments m 
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4. In the opposite case Ti < T* there is a range of temperatures Ti < T < T* where fugacity of single charges Y is 
much smaller than one. This leads to the increase of phase transition temperature, but it still has to remain smaller 
than . The difference with the case of Y ~ 1 is that for Y <C 1 the concentration of free single charges n\ remains 
small even at the temperatures considerably higher than T\ . In the region T\ < T < T2 it can be estimated with the 
use of the standard Debye-Huckel approximation which gives for n\ the self-consistent equation: 
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(where K = T/2Ec) the solution of which for small n\ can be expressed as 
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where a — 0.276.... The main effect of Coulomb interaction is seen in that it produces singularity at T — in the 
exponent in Eq. (|J) . If the shift of the phase transition temperature is taken into acccount Tj should be substituted 
by T\. Note that second and third terms in the numenator of the exponent almost canceal each other in the relevant 
range of parameters. 

For D(T) 3> T the screening of the interaction is noticable only on the large scales and the concentration of free 
double charges also remains small. On the other hand at the temperature T2 = 4T\ « Ec / tt the free double charges 
have to appear even when Y = due to mutual influence of pairs of double charges (cf. with Ref. || ) . That means 
that for D(T) ^> T in the vicinity of T2 there occurs a crossover characterized by a prolifiration of free double charges. 

Close to the transition temperature T\ the self-consistent approximation is no longer valid and more advanced 
methods should be used. It is easy to show that when only the single and double charges are taken into account, 
the Coulomb gas model described by partition function (Q) becomes isomorphic (in continuous approximation) to the 
sine-Gordon model defined by the Hamiltonian: 
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The renormalization group equations for the Hamiltonian (0) can be found in Ref. p7| . As can be expected their 
solution shows that for temperatures lower than T2 the presence of double charges (in form of neutral pairs) does not 
introduce any qualitative changes. In the close vicinity of T\ the temperature dependence of n\ deviates from the 
self-consistent result (0) and follows the standard BKT critical behaviour 101 with 
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where b is of the order of unity. The array's linear dc resistance should be inversely proportional to the density of 
free charges n±. Taking into account Eqs.(^|, ^), we get an estimate for this resistance at the temperatures near T±: 
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where R\ is inversely proportional to the probability of tunnelling event which is only weakly dependent on the 
temperature: R\ ~ R n A4, where R n is the normal-state tunnelling resistance (cf. ]T^|). 

5. The representation (^) is also useful for comparison between our results and Efetov's treatment of screening 
by quasiparticles jl2|], which can be expressed just as the replacement of the original capacitance matrix by the 
"effective" one defined (at T < A) as 
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Let us now formally expand the second term in Eq.(^) up to the second order in 9, and neglect the rest of terms. One 
can easily see that the expression obtained in this way for the effective interaction between 2e charges (generated by 
expansion of the partition function in powers of the cos 29 term) would coinside with the one obtained by inversion 
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of the effective capacitance (|f 0|). Physically the above formal operation would mean neglect of the discrete nature 
of electric charge, which could be reasonable if charge transport between islands could proceed via some "classical" 
channels able to provide charge in continuous amounts (like charging of macroscopic electric capacitor by external 
voltage source). It is not the case for submicron arrays with tunnel junctions where charges of islands can change only 
by le quantum (which is just reflected by periodic nature of the second term in Eq. ([?])). However, basic qualitative 
feature of Efetov's result - screening of Cooper pair charges by normal exitations - rcmanins valid, in spite of the 
absense of any simple notion like an effective capacitance matrix. 

6. It follows from our results, that in order to observe some growth of effective activation energy [defined as 
E a = dh\R(T)/d(l/T)\ with T approaching T 2 from above, one needs to use array with the Coulomb energy Eq few 
times below than A, so that Ti <T*. This conclusion is in complete agreement with recent experimental data JlT| , 
where some moderate growth of E a (T) in the temperature range around 0.2-0.3 K was observed. This experiment 
differs from few previous ones of the same type (where constant E a was observed) by lower values of Ec and a bit 
higher reported A. On the other hand, we do not agree with the interpretation of that E a (T) growth as precursor of 
the BKT transtion at T2, given in [fill for their SC arrays. As follows from our results, no such a transtion exists at T2, 
which agrees with a rather modest (compared to BKT behaviour) growth of E a (T) observed in JllJ at T ~ T2. Note 
that the agreeement between the data reported in [|lTj| for normal arrays and the expected genuine BKT transition at 
Ti is much better than the above-mentioned comparison for SC arrays. 

The above theoretical results point out that any analysis of experimental data on superconductor-insulator transtion 
in artificial arrays of superconducting islands (as well as in the dirty thin films near SC-I transtion) should take into 
account an existence of a characteristic temperature scale of the parity effect T* (note that T* is magnetic-field 
dependent and strongly suppressed by the fields of the order of H c2 ). In particular, the behavior of I-V characteristics 
in the intermediate temperature range T* < T < A cannot be unambigously related to the genuine ground-state 
properties of the system, as is examplified by non-monotonous R(T) behaviour observed in for a nearly critical 
ratio Ej/Ec- We interpret this unusual behaviour as follows: at moderately low T screening by quasiparticles is 
effective and reduces Coulomb repulsion of Cooper pairs, leading to the decrease of R(T) behaviour; at still lower T 
this screening is gone, Coulomb repulsion increases and effective ratio Ej/E^* enters the "insulative" part of the 
phase diagram, leading to the increase of R(T) at further T decrease. 
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